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ABSTRACT 

Electron  Beam  welding  (EBW)  is  modeled  using  a  basic  Kaplan  keyhole  formation  as  in  laser 
beam  welding.  As  high  penetration  depth  to  width  ratios  can  be  obtained  in  this  process,  researchers 
and  engineers  are  working  to  enhance  this  vacuum  based  technology  of  component  joining.  This 
welding  technique  has  several  applications  in  fusion  reactor  fabrication  of  various  components  like 
vacuum  vessel  and  test  blanket  modules  with  different  types  of  steels.  This  study  is  focused  on 
simulating  and  modeling  EBW  in  terms  of  the  heat  flux  flow  with  surface  heat  on  bead  area  in  the 
form  of  Gaussian  distribution.  When  the  beam  forms  keyhole  plasma,  the  heat  flow  conducted  to  the 
base  metal  is  modeled  as  a  frustum  distribution.  Unlike  laser  beam  welding,  in  EBW,  a  portion  of  the 
heat  would  leave  out  to  ambience  through  the  lower  end  of  the  bead,  where  this  amount  of  heat  is  to  be 
apportioned.  Heat  flow  mechanism  is  calculated  using  Peclet  number  from  where  the  empirical 
formulae  power  distribution  is  determined.  Peclet  number  variation  is  calculated  as  a  function  of 
radial  distance,  penetration  depth  and  distance  from  wall.  The  values  obtained  may  not  be  directly 
used  in  analysis  but  a  correlation  is  obtained  which  assists  in  formulating  a  finite  element  analysis. 
The  welding  process  is  evaluated  using  temperature  profile  distributions  in  Weld  Zone  (WZ)  and  Heat 
Affected  Zone  (HAZ)  using  Ansys  software.  SS304  weld  material  is  used  and  temperature  dependent 
thermal  and  structural  properties  are  considered.  The  output  obtained  is  the  thermal  isotherms  around 
the  weld  to  estimate  the  surface  heat  flux  distribution  on  the  surface  and  through  the  keyhole.  The 
scope  of  applications  where  the  heat  flow  distributed  in  %  contribution  would  provide  guidance  in 
structural  joints  as  in  nuclear  reactor  fabrication  applications. 
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1.  INTRODUCTION 

Electron  beam  welding  is  widely  used 
in  technological  applications  where  the 
structural  integrity  is  highly  demanding  in 
terms  of  very  low  distortions  and  low  heat 
affected  zone  with  good  mechanical  properties. 
Electron  beam  technique  employs  high  power 
density  supply  to  the  localized  portion  of  the 
structure  with  fine  control  and  produces  low 
residual  stress  with  narrow  weld  bead 
formation.  Recent  applications  of  the 
structural  components  in  the  nuclear  sector  like 
fusion  and  fission  reactors  have  shown  keen 
demand  in  the  electron  beam  applications  due 
to  the  stringent  working  conditions,  thermo¬ 
mechanical  specifications  and  integrity  of  the 


complex  structures  [1,  2].  In  spite  of  the 
expensive  technique,  the  application  is 
explored  for  demonstrative  purposes  as  well  as 
the  final  engineered  products.  Analysing  the 
state  of  the  residual  stress  and  temperature 
profile  at  several  stages  will  be  very  helpful  to 
design  the  product  features  like  the  influential 
characteristics,  metallurgical  properties  and  the 
liquid  melt  pool  behaviour  during 
solidification  with  the  material  properties.  It  is 
a  challenging  task  to  implement  the  conditions 
such  as  the  high  temperature  properties, 
structural  status  and  the  models.  This  paper  is 
focussed  on  the  analysis  of  the  electron  beam 
weld  process  pertaining  to  SS  304  L  plates 
using  Peclet  number  analysis  where  the 
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empirical  thermal  heat  conduction  to 
convection  ratio  is  performed.  Also  as 
observed  in  the  similar  power  density 
application  like  laser  beam  welding,  the 
differences  are  identified  and  applied  to  the 
analysis,  based  on  the  thermal  isotherms  with 
the  heat  source  input  models  such  as  Gaussian 
beam  source  and  frustum  type  heat  source 
models  [3-5].  The  electron  beam  profile  during 
the  plasma  key  hole  spread  distribution  is 
analysed  with  equivalent  full  Gaussian,  partial 
frustum  and  Gaussian  source  heat  flow  models. 
The  results  are  analysed  in  terms  of  thermal 
distribution  towards  the  weld  and  heat  affected 
zone  planes. 

In  EBW,  free  electrons  of  charge 
1.6X10'19C  are  accelerated  using  electric  field 
and  focused  on  joining  metal  space  as  shown 
in  figure  Al.  Higher  the  electron  beam  weld, 
higher  is  the  metal  penetration  and  material 
melting.  A  vacuum  space  is  created  around  the 
surface  to  avoid  beam  distortion  and 
dissipation  and  to  maintain  focussing.  In  the 
presence  of  air  or  gas,  collision  with  gas 
molecules  would  create  ionization  of 
molecules  and  destroy  cathode.  Vacuum 
prevents  chemical  corrosion  and  eliminates 
blow  holes  and  bubbles  that  might  happen 
even  when  inert  gases  are  used  as  in  LBW.  Just 
like  LBW  the  focused  beam  restricts  fusion 
zone  area  and  hence  has  very  high  power 
density,  up  to  1010W/m2  where  area  is  of  109 
m  .  Initially,  the  beams  transfer  heat  to  the 
surface  of  fusion  zone.  When  heat  gets 
accumulated,  it  is  penetrated  to  few  microns 
depth  and  the  molten  metal  is  exposed  to 
severe  heat  vaporisation  and  gradually  a 
vapour  cover  is  formed  over  the  entrapped 
beams.  Keyhole  formation  occurs  as  how  in 
the  LBW  process  and  there  will  be  more 
efficient  usage  of  beam  energy  to  be  used  in 
the  fusion  process.  In  this  study,  the 
contribution  of  heat  on  the  surface  and  the 
inner  keyhole  is  distributed  and  analysed  by 
the  use  of  Peclet  number.  Different  ratios  of 
surface  to  keyhole  power  distribution  are  noted 
from  the  estimated  power  and  its  effect  on 
temperature,  distortion  and  residual  stresses  are 
determined. 

2.  FORMATION  AND  SIMULATION  OF 
WELDING  PROCESS 

The  heat  energy  equations  are 
referenced  in  many  articles  including  an 
isotropic,  conductive  material  having  equal 


coefficient  of  conductivity,  kx,  ky,  kz(W/mK)  in 
all  three  chosen  orthogonal  co-ordinates. 
Equation  (2.1)  gives  the  heat  energy  in  the 
weld  area  with  temperature  (T(K))  obtained 
both  in  spatial,  x,  y,  z  (m)  and  temporal  (t 

Q  Q 

(sec)),  terms.  Q  (W/m  or  J/m  s)  is  the  net  heat 
from  the  input  and  the  losses  in  the  form  of 

Q 

convection  and  radiation.  Density  (p  (kg/m  )) 
and  specific  heat  capacity  (c),  give  the  right 
hand  terms  of  how  much  heat  is  retained  with 
respect  to  time  in  the  material  and  how  much  is 
taken  away  with  the  velocity  of  welding  (v 
(m/s)).  The  essential  boundary  conditions 
given  throughout  the  body  are  T0(x,  y,  z,  0)  at 
time  zero  or  at  the  starting  of  the  weld.  In 
addition,  the  natural  boundary  conditions  have 
to  be  applied  that  consist  of  normal  conduction 

(Kn  ^)  heat  flux  (q),  convection  (h  (T-T0))  and 

the  radiation  term  (cje(T4  —  Tq  )). 


kx5!Z  +  ky^I  +  kz5!Z  + 


52x 


S2y 


52z 


(2.1) 


Together,  the  boundary  conditions  are 
summed  up  as  in  (2.2): 

Kn  —  q  +  h(T  —  TO)  +  as(T4  -  TO4)  =  0  (2.2) 


When  symmetric  and  insulation 
boundaries  are  considered  as  adiabatic,  with  no 
heat  flowing  through  the  surface,  they  are 
obtained  by  making  convection  zero  and 
conduction  zero  from  the  surface,  where,  Kn  is 
the  thermal  conductivity  normal  to  the  surface 
(W/mK),  h  is  the  convective  heat  transfer 
coefficient  (W/m  K),  8  is  emissivity  of  surface 
radiating  and  a  is  the  Stefan  Boltzmann’s 
constant  (5.67*  10"8  W/m2K4).  When  it  is 
difficult  to  use  radiation  boundary  condition,  it 
is  combined  to  convective  heat  flux  by  using  a 
modified  coefficient  (hY)  for  hot  rolled  steel 
plates  with  an  error  of  about  5%  as  shown  in 
equation  (2.3) 

hr  =  2.4  *  10“3  (E  T161  (2.3) 


Radiation  inclusion  will  increase 
solution  time  by  about  three  times  and  hence  is 
combined  with  convection. 

2.1.  Peclet  number  definition  and 

information 

Peclet  number  gives  thermal  flow  in 
liquid  and  from  liquid  to  solid  state.  Peclet 
number  at  top  surface  of  weld  is  Pe(0)  at  Z=0 
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vertical  distance  and  is  Pe(z)  at  any  depth  z.  If 
the  bead  is  conical  or  of  frustum  shape,  as 
depth  (d)  increases,  diameter  of  bead  width  (D) 
decreases.  Peclet  number  is  a  function  of  weld 
width,  weld  speed  (v),  and  thermal  diffusivity 
(a)  where,  a  is  the  radius  of  heat  source  power 
and  m  is  an  empirical  multiplier  as  shown  in 
equation  (2.4) 

pe(z)  =  pe(0)  (l  —  ^  and  pe(0)  =  v  *  a/2  (2.4) 

Because  of  heat  flow,  the  amount  of 
power  absorbed  (PI)  in  the  inner  side  of  weld  is 
given  as  shown  in  equation  (2.5) 

PI  =  t(ml  *  Kmol)  *  (Tv  -  T0)(2.1995  + 
3.1481Pe(0)  -  0.16647Pe(0)2  + 
0.01152Pe(0)3)  (2.5) 

Here,  t  is  the  thickness  of  plate  to  be 
welded  and  ml  is  empirical  value  to  take  care 
of  variation  in  velocity  and  thermal 
conductivity.  Kmol  is  the  thermal  conductivity 
of  molten  metal,  Tv  is  heat  source  temperature 
and  TO  is  initial  temperature. 

Pt  =  Pc  +  Pi  (2-6) 

where  pt  is  the  total  power  consisting  of 

surface  induced,  pc  and  pj  are  from  the  molt 
metal  laterally  transferred  to  the  solid  base 
plate  along  with  the  thickness  by  conduction 
and  in  the  presence  of  plasma  key  hole  as 
given  in  equation  (2.6). 

In  this  study,  different  combinations  of 
pc  and  pj  are  analyzed  with  the  effects  of 
temperature  change  and  distortion  and  residual 
stresses  that  occur.  The  power  calculated  for 
different  pt  are  shown  in  figure  A10.  The 
variation  of  power  for  different  bead 
temperatures  is  shown  at  different  distances 
from  the  center  of  the  weld. 

2.2.  Finite  element  for  simulation 

The  heat  equations  can  be  represented 
in  tensor  form  so  that  the  elemental  transient 
heat  equation  is  obtained  and  then  summed  to 
get  the  system  equation  which  is  analysed  with 
time  as  in  equation  (2.7) 

LK(T)J{T)  +  [C(T)]{T)  =  (Q(T)}  (2.7) 

where  K  is  a  temperature  dependent 
conductivity  matrix  and  C  is  the  temperature 
dependent  capacitance  matrix  based  on 
specific  heat,  where  their  product  with  rate  of 
temperature  results  in  heat.  The  above  equation 


can  be  solved  numerically,  with  standard  FEM, 
Crank  Nicholson  or  Euler  time  integration 
models.  Assuming  an  initial  temperature  (Ti), 
K,  C  and  Q  are  calculated  at  that  temperature 
and  the  next  temperature  T,  where  i+1  is 
obtained.  Again  K,  C  and  Q  are  calculated  and 
temperature  at  next  temperature  interval  is 
calculated.  The  iteration  is  continued  for 
temperature  convergence  or  heat  flux  values. 
This  is  a  procedure  for  transient  finite  element 
analysis.  In  the  present  study,  the  work  is  done 
using  Ansys. 

2.3.  Finite  element  model 

The  finite  element  model  of 
dimensions  40mm  x  150  mm  x  5mm  is  used. 
AISI  304  austenitic  stainless  steel  material  is 
considered  for  simulations  to  be  carried  out. 
Convection  is  applied  on  all  the  surface  of  the 
plate  except  on  the  heat  applied  area.  In  the 
this  study,  AISI  type  304  stainless  steel  is  used 
as  it  has  many  advantages  such  as  low  thermal 
conductivity,  high  resistance  of  corrosion  and 
high  stability  at  elevated  temperatures.  Thus 
SS304  material  is  widely  used  in  numerous 
industries  like  nuclear  industry,  chemical 
plants,  aeronautical  and  specialized  pipe 
industries. 

The  properties  of  a  typical  stainless 
steel  sheet  and  the  temperature  dependent 
thermal  properties  for  AISI  304  stainless  steel 
material  are  given  in  table  1  and  2  respectively. 


Table  1. Mechanical  properties  of  AISI  304  Steel 


Tensile 

strength 

Yield 

strengt 

h 

Den 

sity 

Melti 

ng 

point 

Thermal 

conductivit 

y 

% 

elonga 

tion 

515 

MPa 

205 

MPa 

800 

0 

kg/ 

m3 

1400- 

1450° 

c 

16.2  W/m° 
K  at  100°c 

20-40 

Table  2. Temperature  dependent  thermal  properties 
for  AISI  304  Austenitic  stainless  steels 


S.No 

Temp 

(K) 

Thermal 
conductivity, 
W/m  °  K 

Density, 

Kg/m3 

Specific 

heat, 

J/KgK 

1 

200 

11 

8200 

350 

2 

400 

15.5 

8000 

400 

3 

600 

19 

7800 

440 

4 

800 

22.5 

7600 

550 

5 

1000 

26 

7500 

590 

6 

1200 

30 

7400 

610 

7 

1400 

34.5 

7350 

640 

8 

1600 

39.5 

7300 

680 

9 

1800 

44 

7200 

720 

10 

2000 

47 

7200 

760 
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3.  THERMAL  ANALYSIS 

Thermal  analysis  has  been  carried  out 
with  constant  heat  flux,  where  the  thermal  load 
is  applied  at  a  time  on  the  weld  area.  Load  is 
applied  for  first  10  sec  and  the  temperature 
attains  fusion  in  40  seconds.  The  time  step  is 
progressively  increased  up  to  time=1000sec  to 
see  the  plate  cooling  down  to  ambient 
temperature.  In  the  present  work,  Finite 
Element  Analysis  (FEA)  of  single-pass  butt¬ 
welding  has  been  carried  out  with  constant 
heat  flux.  A  simple  butt-joint  welding  whose 
welding  parameters  are  consistent  to  those  of 
Friedman’s  model  with  heat  input  Q=2000W  is 
considered  and  has  been  simulated  using 
ANSYS.  The  present  thermal  Ansys  is 
conducted  using  element  type  SOLID70.  The 
element  is  applicable  for  three  dimensional, 
steady-state  or  transient  thermal  analysis.  The 
element  can  also  compensate  for  mass 
transport  heat  flow  from  a  constant  velocity 
field.  In  this  analysis,  element  SOLID70  is 
replaced  by  a  three-dimensional  (3-D) 
structural  element  SOLID45. 

This  element  type  has  a  three- 
dimensional  thermal  conduction  capability  and 
eight  nodes  with  single  degree  freedom 
(temperature)  at  each  node.  The  element  is 
defined  by  eight  nodes  having  three  degrees  of 
freedom  at  each  node  (translations  in  the  nodal 
x,  y  and  z  directions).  The  element  has 
plasticity,  creep,  swelling,  stress  stiffening, 
large  deflection  and  large  strain  capabilities  as 
depicted  by  ANSYS.  The  geometry  and 
meshed  model  with  tetrahedral  shape  with  a 
volume  mesh  of  size  0.02  were  shown  in  figure 
A2.  Figure  A3  and  A4  show  the  geometry  (3mm 
thick)  and  mesh  model  used  for  analysis 
respectively. 

[6-8]  provides  expressions  and 
numerical  models  for  different  heat  fluxes 
including  Gaussian,  where  heat  flux  is  in  the 
shape  of  a  normal  distribution  used  for  surface 
heat  on  fusion  zone.  Also,  frustum  is  a  part 
obtained  by  removing  a  part  of  the  apex 
portion  from  a  cone,  where  this  distribution  is 
used  for  keyhole  heat.  A  combination  of 
Gaussian  on  the  surface  and  a  conical  heat 
distribution  in  the  key  hole  is  recommended 
for  laser  welding  process  as  the  bead  ends  in 
the  thickness  of  the  material.  It  was  appropriate 
to  use  a  conical  heat  flux  distribution  in  the 
lateral  surface  of  material.  In  case  of  EBW, 
though  the  heat  flux  distribution  is  same  as 
EBW,  the  power  of  the  beam  makes  the  beam 


open  at  the  bottom  of  weld.  This  is  better 
modeled  by  a  frustum  distribution  of  heat  flux 
in  the  bead  area  or  penetration  area  of  the 
thickness  of  the  plate  instead  of  a  cone.  Here, 
the  combination  of  Gaussian  and  frustum  are 
used  for  heat  flux  model. 

4.  RESULTS  AND  DISCUSSION 

Thermal  results  with  heat  flux  values 
are  obtained  from  Peclet  number  and  used  in 
EBW  heat  flux  models.  Mechanical  properties 
of  AISI  304  steel  are  shown  in  table  1.  As 
temperature  has  an  effect  on  both  structural 
and  thermal  properties,  its  variation  is 
incorporated  in  the  analysis  as  shown  in  table 
2.  The  variation  of  mechanical  and  thermal 
properties  is  incorporated  in  both  thermal  and 
structural  analysis  of  the  sequentially  coupled 
procedure  used. 

The  temperature  distribution  was 
found  to  be  linear,  as  represented  in  figure  A5. 
X  axis  is  a  ratio  of  surface  to  lateral  heat  and 
since  the  lateral  heat  is  encompassed  in  the 
keyhole,  convection  in  molten  metal  helps  in 
fusion  welding  and  excess  heat  is  transferred  to 
the  wall  and  conducts  freely. 

The  temperature  distribution  profile 
with  color  coding  is  shown  in  figure  A6.  For  a 
50:50  equal  distribution  of  surface  and  lateral 
heat,  the  temperature  at  the  base  and  the  bead 
is  about  470°K  and  2586°K  correspondingly. 
Similarly,  in  the  40:60  distribution  of  heat  ratio 
of  surface  lateral  heat,  it  is  found  to  be  492°K 
and  276 1°K  respectively.  In  figure  A7,  the 
temperature  change  at  bead  variation  with  time 
is  shown.  It  is  shown  for  30:70  ratio  of  surface 
to  lateral  heat,  where  the  maximum 
temperature  is  2800°K  for  about  40  seconds 
and  after  1000  seconds  it  is  850°K. 

4.1.  Structural  results 

The  structural  results  are  shown  for  the 
power  (30%  Gaussian,  70%  Frustum) 
distribution  of  the  model.  Due  to  the  variance 
in  the  temperature  gradient,  the  material 
properties  are  given  in  the  model  for  elevated 
temperatures.  A  stress  acting  normal  to  the 
direction  of  weld  bead  is  known  as  a  transverse 
residual  stress.  Distortion  is  shown  in  figure 
A8,  where  the  minimum  distortion  is  at  the 
base  metal  and  the  maximum  distortion  is  at 
the  weld  bead  is  0.2mm.  The  minimum  V  on 
Mises  stress  are  at  the  base  plate  and  the 
maximum  Von  Mises  stress  occurs  at  clamping 
points  305MPa,  otherwise  it  is  within  llOMPa 
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at  all  other  points  as  shown  in  the  figure  A9. 
The  temperature  near  the  weld  bead  and  heat 
affected  zone  rapidly  changes  with  distance 
from  the  heat  source.  This  shows  more  stress 
value  in  the  weld  bead  area  and  gradually 
decreases  from  center  line  to  the  base  plate 
end.  Both  distortion  and  stresses  are  in 
acceptable  limits  in  EBW. 

5.  CONCLUSIONS 

An  electron  beam  welding  model 
using  the  keyhole  formation  and  the  Peclet 
number  is  discussed. 

•  The  heat  distribution  is  divided  among 
the  two:  surface  using  Gaussian  and 
Frustum  distribution  for  lateral  heat. 

•  A  variation  in  the  heat  distribution  is 
tried  in  which  the  surface  to  lateral 
heat  equals  from  50:50  to  30:70. 

•  A  linearly  increasing  distribution  with 
temperature  getting  better  at  bead  with 
increase  in  lateral  heat  entrapped  in 
keyhole  is  beneficial  to  obtain  weld. 

•  Weld  temperature  would  be  obtained 
in  about  40  seconds  and  rest  of  1000 
seconds  would  be  for  cooling  down  to 
ambient  temperature. 

•  Distortions  and  stresses  are  stable  with 
heat  flux  distribution  changes  and  are 
well  within  the  acceptable  strength 
limits  at  205MPa. 
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Figure  Al.An  electron  welding  beam  generation  equipment. 
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Figure  A2.Mesh  model  used  for  analysis  and  element  model  used 
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Figure  A3. Geometry  of  the  model  (3mm  thick) 


Figure  A4.Mesh  model  used  for  analysis 
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Figure  A5. Power  vs.  temperature  distribution 
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Figure  A6. Temperature  distribution  of  50%  Gaussian  -50%  Frustum 


Figure  A7a. Temperature  distribution  (40%  Gaussian-60%  Frustum  of  power) 
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Figure  A7b. Temperature  distribution  of  30%  Gaussian  -70% 

Frustum 
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Figure  A8. Distortion  for  30%  Gaussian  -  70%  Frustum 


Figure  A9.Von  Mises  stress  for  30%Gaussian  -70%  frustum 


Figure  At 0. Power  estimation  using  Peclet  number 
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